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TECHNICAL MEMORANDUM 119T 


INVESTIGATIONS ON THE STABILITY, OSCILLATION, AND STRESS 
CONDITIONS OE AIRPLANES WITH TAB CONTROL 
FIRST PARTIAL REPORT - DERIVATION OF THE EQUATIONS OF 
MOTION AND THEIR GENERAL SOLUTIONS* 

By B. Filzek 

For the design and the construction of airplanes the control Is 
of special significance, not only with regard to the flight— 
mechanical properties "but also for the proportional arrangement 
of wing unit, fuselage, and tail unit. Whereas these problems 
may he regarded as solved for direct control of airplane 
motions, that is, for Immediate operation of the control surfaces, 
they are not clarified as to oscillations, stability, and stress 
phenomena occurring in flight motions with indirect control, as 
realized for instance in tab control. Its modus operand! 1 b 
based on the activation of a tab hinged to the trailing edge c£ 
the main control surface. Due to lift and drag variations, 
moments originate about the axis of rotation of the main control 
surface which cause an up— or— down floating of the main control 
surface and thus a change in the direction of the airplane. 

Since this tab control means flying with "free control surface", 
the treatment of this problem should provide the basis on which 
to Judge stability, oscillation, and stress data. 

Aside from the discussion of free oscillations and forced 
oscillating conditions, this report points out methods to deal, 
above all, with starting conditions. The required mathematical 
expedients are suitable, not only for the solution of the given 
problem, but also of similar problems, and are therefore treated 
In general formulation. 

The present report Is to represent a contribution toward the 
clarification of the problems arising and, first of all, to treat 
the longitudinal motion of an airplane. 

*"Untersuchungen uber die Stabilitats— , Schwingungs— und Beanspruchungs— 
verbaltnisse von Flugzeugen mit Hilf srudersteuerung. 1. Teilbericht. 
Herleitung der Bewegungsgleichungen und ihre allgemeinen Losungen." 

V Zentrale fur wiBsenschaftliches Berichtswesen der Luftfahrtforschung (ZWB) 

Forschungsbericht Nr. 2000, October 24, 1944. This translation is the 3 
first partial report (l. Teilbericht) of an investigation made up of two 
„ parts, the second part, FB 2000/2 (2. Teilbericht), of which is NACA TM II 98 . 
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I. INTRODUCTION 
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The control of airplanes up to flying weights of about 20 tons is 
possible with direct control, that Is, direct operation of the control 
surfaces by the pilot. However, for the modern large aircraft of about 30 
tons flying weight and still more for the large aircraft to be expected 
in the future, the power required for direct activation of the control 
surfaces Increases so greatly that human muscle power is no longer sufficient. 
One is obliged to use indirect controls, the necessary power being supplied 
for Instance by a servomotor. The tendency will be to manage with a 
minimum of power. 

The present report Is to represent a contribution toward clarification 
of the problems arising and, above all, to treat the longitudinal motion 
of an airplane . 


II. SYMBOLS 

i 


The symbols used correspond to the German standards DIN L 100 second 
edition, July 1939 } moreover, the following symbols are chosen (cf. also 
fig. 1): 
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m 

m L 

Mg 


Sg 

S 3 

®1 


e. 




^Hr 

F Hh 

^h 


Ati 


bmax 


Let 

l H1i 

C MF 

c nH 


mass of the entire airplane 

partial mass fuselage without control surfaces and tabs 
partial mass control surfaces without tabs 
partial mass of tabs 

center of gravity of the whole system. 

center of gravity of the partial masses m^ 

center of gravity of the partial masses mg 

center of gravity of the partial masses m^ 

distance of the center of gravity of the partial 

mass m][_ from. the center of gravity of the whole system 

distance of the center of gravity of the partial 

mass Mg from the center of rotation of the main control 
surface 

distance of the center of gravity of the partial 
mass from the center of rotation of the tab 

moment of inertia of the entire airplane about y-axis 


r 

F S1 

moment of inertia 
through S-|_ 

of the partial mass 

m l 

about y-axis 


% 

moment of 

inertia 

of the partial masB 

mg 

about 

s 2 

-4 

% 

moment of 

inertia 

of the partial mass 

m 3 

about 

S 3 


moment of Inertia of the elevator about axis of rotation 
of the main control surface 

moment of inertia of the tab about axis of rotation of tab 
deflection of the tab 

deflection of the tab required for attaining the safe 
multiple of load 

chord of the elevator + tab 

chord of the tab 

coefficient of the air— force moment of the airplane without 
horizontal tail referred to center of gravity of the 
airplane 

coefficient of the normal force of the horizontal tail 
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coefficient of the horizontal tail moment referred to the 
axis through the aerodynamic center of the horizontal 
tail parallel to the y-axis 

C rH coefficient of the elevator moment referred to the axis 

of rotation of the main control surface 

C-rTT>i coefficient of the tat moment referred to the axis of 

rotation of the tah 

angle of attack at the location of the horizontal tail 

% mean dynamic pressure at the location of the horizontal tail 

surface 


mean velocity at the location of the horizontal tail surface 

mean dowrrwash angle at the location of the horizontal tail 

ratio of the air force damping of the entire airplane with 
respect to lateral axis to the air force damping of the 
horizontal tail 

If similar symohls appear in the individual paragraphs (for 
instance, Aa = angle-of-attack increment, AM^ = moment increment, A = dowcv- 
wash angle, A k = determinant. A* = South's discriminant), they are dis- 
tinguished hy indices and their significance follows from the respective 
connection. 


V H 

A 


III, PRESUPPOSITIONS 


In order to circumscribe the range of validity of this report, the 
necessary presuppositions shall he mentioned in advance; at the respective 
places they will he especially emphasized .- 

1. Unsteady air force influences are not taken into account. 

2. Flight motions are possible only in the XgZ g -plane. 

3 . Method of smal 1 oscillations. 

4. Fully compensated control surfaces, that is, the centers of 

gravity lie on the axes of rotation of the control surfaces. 

5* Omission of the mass couplings and hprizontal tail surface 
forces in the force equations. 
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6 . Omission of the mass couplings and of the force couple at the 

horizontal tail at C aH =- 0 (zero moment ) in the equation of 

moments about the lateral axis. 

7. The velocity along the flight path v for the phenomenon under 

consideration is assumed to he constant. 

8 . Air forces and mome nts are assumed to he linear functions of their 

variables . 

9 . The steady state about which the system oscillates is assumed to 

be horizontal flight. 

in . The dynamic pressure ratio q^/q. ani ^ downwash factor c&/cta. 
are regarded as constant during the flight motions. 


IV. DERIVATION OF THE EQUATIONS OF MOTION 


For obtaining the equations of motion Lagrange's method will be 
used which permits the required equations of motion of a system 
with n degrees of freedom to be found from the following generally 
valid relation: 


d_ 

dt 


St 



(1) 


Since thiB report deals with accelerated motions, it is useful to select 
for the determination of the kinetic energy T a nonaccelerated system 
of reference. For the present problem the coordinate system fixed to the 
ground may be regarded as such a system of reference. In 'that system, 
in its most general formulation, one would have to consider for an air- 
plane in free flight the degrees of freedom 

1 . translation of the center of gravity of the airplane in 

X— direction, 

S 

2. translation of the center of gravity of the airplane in Y -direction, 

O 

* 

3 . translation of the center of gravity of the airplane in 

Z —direction, 
g 

4 . rotation of the longitudinal inclination of the airplane about 

the center of gravity of the airplane d, 

5 . rotation of the lateral inclination of the airplane about the 

center of gravity of the airplane <p , 
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6. rotation of the compass angle about the center of gravity of 

the airplane y, 

7 . rotation of the elevator about axis of rotation of the main 

control surface tj. 

However, far the consideration of the longitudinal motion of an 
airplane due to symmetric elevator operation one may assume that the 
airplane executes motions only in the XgZg— plane. Thus no motions 

take place in the coordinates Yg, cp, and \jr and the system is reduced 
to the degrees of freedom g > iS, and T) • 

For the system sketched in figure 1 the kinetic energy assumes the 

form 



The velocity coordinates vg r of the partial masses m^. result from the 
position coordinates X^, by differentiation with respect to time 

as the relations 


J 
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r s 2 = X g 2 + Z g 2 


Tg ^ 2 = Xg 2 + Zg 2 + S ^ 2 — 2 ©^ XgSin a + Z g cos a 


T s 0 2 = * g 2 + 2-g 2 + r H 2 ^ 2 + ®2 2 ^ + ^i) 2 + 


XgSin -a + ZgCos aj 


+ 2e 2 (a + ti) |^g sin (iS + *]) + Zgcos (,g + ri)J 
+ 2r2©2^(^ + T[) jsin a ein (a + tj) + cos a cos (a + T})j 

^S 3 2 ■ Xg 2 + Zg 2 + rg2.a 2 + ^Lpr - Lnh^ 2 (a + q) 2 + e 3 2 ($ + n + %)' 

+ Sr^ajXgSin a + Z g cos aj 

+ 2^. - L^OS + rj) [^gSln (a +tj) + Z g cos (a + tj)[ 

+ 2e 3 (^ + ti + t)k) X g Bin(a + ij + + Z g coe^ + i\ + 

+ 2r H (LHr - + n) COB n + 2r H e 3*(^ + i + 008 (ti + ^ 

+ 283^2^. - •S + Tl)(^ + T) + % N ) COS TJk 


X3) 
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If the derivatives are formed according to the rule given in equation (I), 
one recognizes that the equations of motion are no longer linear and thus 
cannot "be treated in a closed mathematical form. Therefore the expres- 
sions are already linearized at this point, that is, one puts cos X = 1 
and sin X = 0 . 

These deliberations together with the summarizing of certain expres- 
sions yield for (2) the form 


T - |{ m ( i g 2 + V 2 ) + V 2 + J Hrn 2 + *Wlh 2 

+ 2*| [F^ + *3*^(1^ - + mgegrg + a^e^] 

+ 2 ^hfHh + m 3 6 3 { r H + ( L Hp - + 2 * Z g ^(“2 + m 3 ) ~ e l^ 

+ 2qii h ji’g h + m^eg^Lgp - LggjJ + 2(d + Tl)^g |™3 ( L Hr “ ^h) + m 2®2 
+ 2111363 + q + T]g) Zgj. 


/(2a) 


If one introduces for the "generalized forces" Q„ convenient symbols 


and the derivatives 


'cff 


- 

c*. 


into (l), the desired equations of motion 


read, if is regarded not as free coordinate, but as disturbance 

function 

Xgm = -A sin 7 — W cos 7 

Zgm + d J^mg + m3) rg + m3 ^Hr — Lm^) — m 1 e 1 + 21.363 

+ q jm.3 (Lgj, — Lgg) + n^egt m 3 e 3j = G- + W sin 7 — A cos 7 

*g [K + m 3) r H + m 3 ( L Hr - L Hh)“ m l 9 l + m 2 e 2 + ^ + * F y 

+ nj^Er + m 3 r H(V- “ ^ + + m 3 e 3 r H] = 

^g m 3 ( L Hr ” %h) + m 2 e 2 + m 3 s 3 J + ^^Hj, + m 3 r H ( L Hr “ L Hh) + m 2 e 2 r H 

+ ^393^ + nFar = M Hr (° 4 [j W* %) 
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One sees that this system contains the eight variables X g J Z g» 7> 
T), a, ag-j and t^. Since there exist here four relations connecting 

the variables, which are given by 

i 3 = a + 7 "1 


and 


ctg = a H (a, d,d ) 
tan 7 = 


> 


( 5 ) 


% = 1 lh( t ) 


(disturbance function) 


the four equations of motion derived in ( 4 ) are necessary and sufficient 
for a complete description of the problem. 


The terms %.e r and rueyXBy, respectively, contained in this general 
formulation of the equations of motion permit — taking the weight moments 
appearing in and Mj^ into consideration — the couplings caused by 

unmass— balanced control surfaces to be included. 


An assumption of fully mass— balanced control surfaces, that is, 

®1 = e 2 = e 3 = 0* does not only signify an essential simplification from 

the mathematical point of view; it is, with respect to eg and e^, 

necessary to ensure freedom from flutter. 


If one considers, moreover, the orders of magnitude of the individual 
sum terms in the corresponding equations, one finds that for rigid 
aircraft part of the coupling terms is of subordinate significance; if they 
are neglected, one obtains the following simplification of the equation 
system: 


Xgm = — A sin 7— ¥ cos 7 
Zgm = G + ¥ sin 7 — A cos 7 

**y = ^ 

’•a’^Hr + mgra^Lar - Lnhj + il^Hr = M Hr 
or, written in the manner customary for flight mechanics. 


(6a) 
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mjl = -W — G sin 7 


mv^2L = A — G cos 7 
at ' 


V = ' m l 


^ ^Hr + m 3 r H( L Hr " L Hh) + *^ F Hr = 

L 1 


As shown by flight measurements and calculations, the first equation of 
( 6 b) affects the variation of the coordinates 7 and •a only Blightly, 
due to the small change ih speed v along the flight path at starting 
conditions. If one disregards, therefore, this equation for the solution, 
the remaining equations with v = const. = v Q are sufficient for the 

description of the problem and one obtains the equations of motion 

mv og- = A — G cos 7 
dt 

V = m l [ (6c) 


F Hr + m 3 r H( L Hr “ -%h^ ^ + ^HrH = 


the further mathematical investigations will be based on them in this shape. 

The first equation of the equations of motion of ( 6 c) deals with the 
forces perpendicular to the plane of the trajectory. 

The second equation represents the condition for moment equilibrium 
about the center of gravity of the airplane, whereas the third equation 
expresses the condition for moment equilibrium about the axis of rotation 
of the main control surfaces. 

For the "generalized forces" Q rj designated in ( 6 c) by A, 
and one makes the assumption that they are linear functions of the 

variables. Then the expressions 

A = ^%F a (7) 


^Bibliography (4) 
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are valid for ( 6 c, l) and ( 6 c, 2). A and M L are, therefore, relatively 
independent of the method of control "by tats. 

In expression of the external moment M^. in ( 6 c, 3), however, the 
form of the tab control obtains decisive importance. 

In figures 2 and 3 the most frequently used tab controls are repre- 
sented schematically. 

The 'direct control by tabs is constructed so that no direct mechanical 
control possibility exists between control stick and tab. Instead, a 
deflection of the control stick actuates a servo-unit fixed to the main 
control surface and thus effects a displacement with respect to the main 
control surface. 

For indirect tab control a deflection of the control stick actuates, 
by a push-pull control, a lever pivoting on the axis of rotation of the 
main control surface. Due to the relative motion of the lever with regard 
to the main control surface deflections of the tabs are effected. If one 
inserts, moreover, a spring between lever and main control surface, the 
result represents a combination of a servotab control and a direct control . 
If one makes the spring infinitely rigid, the system is transformed into a 
direct control by main control surfaces. 

For the control by tabs visualized in figure 2 the expression 
for reads 

w ^°rH ^rH 1 F T % UK _ /«\ 

M Hr + -5^7 + " A l^ (9) 

the factor U)^ is -to be interpreted as damping factor of the main 
control surface. 

If one wantB, on the other hand, to give the expression for an indirec 
tab control (fig. 3 ), the relation 


= ^CrH + ^Fed + “sth + “d 
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must 


M. 


CrH 


Mi 


Ted 


M, 


Sth 


be taken into consideration, with 

= airforce moment of the not moving elevator 

= moment of the spring 
= moment of the tab rod force 


M D 


9 

7 

% 


= moment of the main control surface damping = — ^A^r) 
If in figure 4 

signifies the deflection of the driving lever a 
the deflection of the main control surface 
the deflection of the tab 


i the mechanical gear ratio between tab deflection and 

driving lever a 

there exists between the three angles, due to mechanical coupling, the 
relation 


-ti h = (9 - Tl)i(cp,Ti) (10) 

Due to the deflection of the spring caused by the activation of the 
lever a, a moment of the following magnitude arises 

M Fed = C Fed f ( c Pj 1 l) ( 1J -) 


with Cped. signifying the moment of the spring about the axis of rotation 
of the main control surface referred to 1°. 

The tab push— rod force Pgth produces a moment about the axis of 
rotation of the main control surface 

* 

M Sth = f 12 ) 

with 


M- 


[Sc 


rHh. 


5c 


Hh 


Sa-H 


a s + 


rHh 


3t] 


+ 





(12a) 
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For small angles the gear ratio may "be put as 


i(cp,T)) = const. = i (10a) 

and 

= (<P - i\) (11a) 

Only variations with respect to a condition of steady motion are 
decisive for the course of the motion. Thus one will introduce for the 
variables the expression 

7 = 7 0 + Ay 

a = a Q + Aa etc . 

One then obtains for instance for the first equation of ( 6 c,l) 
dv ^a 

mY °dt = ^“lo^o + “ G COB <?o + 


and by subtraction of the steady starting condition 

0 = - G cos 7 

aa 


the relation 


, 3c 

mv<y=2i. = — S^oFAa + G sin 7 Q A 7 
dt da 


(13) 


If one treats the expressions for and in a similar manner and 

takes into account that the coefficients c ^tt . CpH may be 

expressed as linear functions of the corresponding angles^ for instance for 



bibliography ( 5 ) 
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and 



orj respectively, talcing (10) to (12) into consideration. 


"jin 


rH + ^cCrEhl 


Sa H 




^rH SCrHh " 

\ St] + br\ j 


+ il AS + 

3 iiJ 


Qtr 

~ C Fed 


rH + 1 ^ c 'rHh 


3C 


Sri, 


^H L H q -^ ~ C Fedj' ACp “ A lr^ 


A Y (915) 


( 1 ), 


the factor ' of the last term of the equation (9h) is to he regarded 

as damping parameter of the main control surface. 


For the angle of attack cw at the location of the horizontal tail 
the expression customary in flight mechanics will he used. 



( 15 ) 


If in (13) the horizontal flight is regarded as starting condition, one 
obtains 



and from the relations connecting the variables (5) result — taking into 
consideration the expression 

d = dg + Ad 


and 


d = 


dd 

dt 


a (Ad) 
dt 


= Ad 


the relations 
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Afl = Aa + Ar 
i = a + 7 = a + C-jjAa 


> 


(16) 


‘■a = a + 7= a+ Cjja j 

With the according summarizations and transformations the relations 
(5) to (16) yield, inserted into (6c), In the conventional form the 
differential— equation system with the two degrees of freedom a and -q 


il + + ^A-^At) + ^ A iot a + ^ A la°- + ^ A laA a = c l AS ( t ) 


> 


(17) 


^A2rjAn + a + ^A2afl- + ^A2aA a = C^^t) 


one has to put for the variation of the disturbance function AS(t) in 
the case of direct tah operation 

AS(t) =An h (t) ■ . 

and In the case of indirect tah operation 

AS(t) = Acp(t) 


The differential equations thus obtained may now 3erve, for a 
prescribed dynamic pressure, for the determination of thd course of the 
movement. Generally, however, it will be desired to treat the entire 
dynamic pressure range in a closed form. If one chooses, therefore. 


t = tf4t as new independent variable and puts At) = At}(t); 

Aa = Ax(t)j Atj^ = At^t) and, respectively, Acp = A<p(t) one obtains 


r\ - rT<L 


and, respectively, a = a’ \[q 


a = a" q 




(13) 


These relations (18) substituted Into (17) permit, after division by 1, 
a representation of the differential equations no longer explicitly 
dependent on the dynamic pressure which may then be written 
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V* + ^AitjV + 


Ai^n + ^Aiop-" + ''■^Aiaa* + ^AiaAa = CiAS(t) 
^°^A2t)Aii 4- a," + ^■ L ^A2a a '* + ^A2aAa = C2 &S(t) 


( 1 ) 


( 0 ), 


> (17a) 


If one takes into consideration that the moment due to in (8) and 

(8a), respectively, remains Hmall in comparison with the other quantities, 
the coefficients may he written 

a) for direct tah control 

^A-j_^ = ^A-^ parameter of the main— control— surface damping 


( 0 ) 


A 1t) = 


d°rH F H L H % 
Sq Fg*. q 


( 2 ) 


Ala 


= ± + ^H^Hr ~ ^h) 


Far 


(!). \{2)^ SlL & g aC rH FgLgrg q 5 

a 8a G ■ 3ag F^ q 


Ala = 


( 0 ) 


Ala = — 


dC’rH ^5% 

3ag F 


Hr 


|i _ & + Sgfi 

8a G 2 3a 


. (% ♦ BjVi 
*■ ^ 


r-^E 


51 

q 


(0) v 

A 2 n - -5T "V r 


( 1 ) A _ ^nH /* 

A 2a - + 


( 0 ) 


A2a = 


5£nH j"l _ c& 
Sag L 3a 


Sc. 




3c 




H°H 


c X = - 


3ag 17 3a T 3a G® 2" “• “■ V qsl “ “q 
3c^g FgLg q 2 



F H r H7^ 


3C 


'MF- 


3a 


-FL 


1_ 

F„ 


S 7h ^ q 


n _ ^°nH F H r H iH 

Wo — 1 . — — 


3rj 


i h r i 




J 


(1Tb) 
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1) for indirect tat control 



In the scope of this investigation it does not signify an essential 
limitation if one regards "both dynamic pressure ratio q. H /q. and dcwnwash 

factor c&/da as constant not only during the flight motion hut also far 
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the dynamic pressure range to "be investigated, so that vith this 
presupposition a linear, inhomogeneous or complete system of differential 
equations with two degrees of freedom is presented for discussion. 


Y. GENERAL SOLUTION' OF LINEAR DIFFERENTIAL 
EQUATION SYSTEMS (DES) 


Knowledge of the solutions of inhomogeneous differential equations 
or, respectively, differential— equation systems is a presupposition for 
the treatment of starting phenomena as they are of importance In flight 
mechanics3 and other technical fields. These solutions are given in the 
most general form for a linear DES with constant coefficients and n degrees 
of freedom with arbitrary disturbance functions and starting conditions. 

1. Solution for Arbitrary Disturbance Functions and Starting Conditions 

One assumes as given the DES with n degrees of freedom 


K=1 p=0 


(19) 


Therein Is: 


K 

= 1, 

2, 

3 ••• 

... n 

the designation of the dependent variables 

V) 

= 1, 

2, 

3 ••• 

• • • n 

the notation of the independent differential 
equations 

P 

= o. 

1, 

2 ... 

. . . m 

the order of derivative of the respective 
dependent variables 


X K ^ = the pr-th derivative of the K — th variable^ is to 

be = Xfl 

(^ a \ DK = the coefficient of the DES which stands for the p— th 
derivative of the k— th variable in the n— th line. 


^Cf. bibliography, (4) and (6) 
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The desired solution of this problem can "be given in a particularly clear 
form by using the Laplace transformation^. The required auxiliary theorems 
are introduced here without any proof. 

If one applies to (19) the Laplace transformation 

F(T)c - P T dT (20) 



and bears in mind that there exists the relation 

cf^lpM(t)|- = P^cjF(t)j- F^^(0)p^ -A-;L 


( 21 ) 


(with the sum to be put equal to zero for p. = 0), (19) is transformed 
into 



or 





with 



+ ^ s oic (p) 


(22) 





(lO* 


, OK- i 


functions at most (23) 
of nt-th degree in p 


g (p) x (^(o)?^- 1 

A—pr uk k K 
|i=l A=0 


that is, functions (2h) 
at most (m — 1)— th 
degree in p 


^Bibliography (8) 
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Therewith the linear DES with constant coefficients is reduced to an 
algebraic problem with the n unknown ^.jx^. • Tile coefficients of the 

unknown 




are rational integral functions at most of the m-th degree 
in p. The absolute terms on the right side of (22) consist of the 
Laplace transforms of the disturbance functions of the DES (19) and 
rational integral functions of at most (m — 1)— th degree in p which 


contain all starting conditions 
of the problem. 


X 


(A) 


(0) with A = 0, 1, 2 (m — 1) 


If one applies to (22) the theorems of the theory of determinants, 
the solution of the unknown oO^X^j- may in the most general form be written 

as follows: 



(25) 


therein D is the coefficient determinant, that is D = 0, the character- 
istic equation of the DES (19) 


D(p) 


f UK (p) 


(26) 


D K in each case the determinant which results from the coefficient 

determinant D if its k— th column is replaced by the sequence of the 
absolute terms of (22) 


e^fp) 


(27) 


Since, however, the determinant D R may be written as sum of two 

determinants one of which contains in the respective column the Laplace 
transforms of the disturbance functions F^(t), the second the sum 

n 

y g UK (p), it will be useful to agree upon the following: 

K=1 

If in (22) all F^(t) = 0 and at least one X K ^^(0) ^ 0 one is 
dealing with a homogeneous component solution; if at least one E (t) 4- 0 

(a) U 

and all Xk' a '( 0) = 0, with an inhomogeneous one. 


With this stipulation (25) may be written also 
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t 


9 



_ ^Klnh. , ^ it ham. 
D D 


(28) 


If one now Imagines the determinants written down and calculated., there is 
1 . D a polynomial of at most (mn)— th degree in p 


mn 


2. D 


Kinh 
since the 


D(p) = (26a) 

a sum of polynomials of at most m(n — l)— th degree in p. 


4-th column contains the sequence of the 


Laplace — transform 


/•{yt)} 


developed with respect to the 


If this determinant i3 

it— th column, one can also write 


P l rfTlh (P) 



with A dk "being the algebraic complement of the element 
polynomial of at most m(n — l)— th degree in p. 



3. D K>ir<tTi is a polynomial of at most m(n — 1)— th degree in p 


Kham 


(p) = 


v>=l 


^ SqkCp) "A ulc (p) 1 (^C k p A (30) 

It=l A =0 


so that after these deliberations (28) may be written in the 
form 


oC^k} 

with A /D and D^^/P representing rational functions in 
the form of real fractions in p. 



J Kham 


(p) 


D(p) 


(28a) 


If the characteristic equation (26a) has only simple roots designated 
by p^(p = 0, 1, 2 (mn — 1)), the rational functions in the form of 

real fractions of the right side of (28) may be represented with the aid 
of decomposition to partial fractions in the following manner 
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W p) . = | s=1 VW k 


d(p) "ir=o“ d ’(p^) p - ?n 


(31) 


and 


^horntp) D tc hqn(PM') 1_ 

p(p) D * (Pll) P “ Pl 


therein 


D'(Pu) 


dD 

dp 


P=Pti 


mn— 1 

J 


^Hpu - Po) 


(32) 


(33) 


with the raised w (^)" signifying that the factor is to be put equal 
to "1” for which u * u- 


If one now introduces the two functions of time 


and 


r (tl -Put 

”« c ’ ' hr 


mn— 1 n , ( p,,) 

(t) e p„t 

^ D’(Pn) 


(3*0 


(35) 


their Laplace transforms will be exactly identical to the decompositions 
to partial fractions represented in (31) and (32). 

Therewith (28) and (28a), respectively, are transformed into 


using the convolution CFaltungsJ theorem 
cCj] F(tjj^jG(t)j F(T)G(t - r)dri 


on fi obtains 



NACA TM 1197 


23 



- T )dr + H^t) 


If one now cancels again the Laplace transformation and takes (34) and 
( 35 ) into consideration, there results as solution of ( 19 ) 


x K (t> = 


_OJS_ 


(Pu) 


■0=1 p=0 Dt (5^) J 

UO 


f mn — 1 D 

M F D (T)e- J Vdr + l KhQm 


(Pti) 


tl=0" D'(p^) 




(36) 


Thus, in the most general form, the solution of a linear DES with 
constant coefficients for arbitrary starting conditions and integrable 
disturbance functions is indicated in closed form. 

It has to be noted that ( 36 ) was derived only under the assumption 
of simple roots p^ of the characteristic equation D(p) =0 of the DES. 


If the solution for multiple roots of the characteristic 
equation D(p) = 0 also is of interest, no particular difficulty is 
encountered if the corresponding expressions of the related decompositions 
to partial fractions are taken into consideration and introduced. 

2. Solution for Constant Disturbance Functions and Arbitrary Starting 
Conditions 

For the practical engineer, however, thiB general, solution gains full 
significance only when he Is given ways, and means to represent the desired 
functions T K numerically for his statement of the problem. The 

solutions may be regarded as known when one succeeds in representing, 

for analytically 'prescribed disturbance functions F^(t), the quadrature 

of the time integral F^(«r)e — ^ VdT In closed form. 

In all cases where a representation In closed form is difficult or, 
respectively, the F^t) are not prescribed analytically but are, for 

instance, known from tests as diagrams, there exists fun d a me ntally the 
possibility of replacing the function Fy(t) with a step function 

according to figure 2 In such a manner that the areas enclosed by the 
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l 

curve Fy(t) or, respectively, by the step function and the abscissa 
are equal. 

According to the mean value theorem of integral calculus, in an 
arbitrary time interval At 



*\j(t)<1t = AtJ^Ctj) = AtF u (t^ + dAt) 


0 < d < 1 


If one performs the limiting process with At— 0 , the step function tends 
toward the function F^(t) according to 



- W) 


that is, one can approximate an arbitrarily prescribed integrable dis- 
turbance function F^t) by according selection of the At with arbitrary 

accuracy. 

One can therefore visualize any arbitrary disturbance phenomenon as 
built up from constant disturbance functions which commence to act 

at the times t^, so that one obtains a very serviceable approximation for 
the component solution X Kinii of the inhomogeneous DES by superposition 

of the partial solutions due to the constant disturbances C^. 

For investigation of starting conditions the explicit form of (36) 
shall thus be given only for 

F 0 (t) = Const. = 

The inhomogeneous part of (36) then obtains the form 
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X 4 V 
Kinh 


n an- 1 

UTZ 

b=l Jls?0 


„ - 1 
u E'(P|J Pp. 


which can still he considerably simplified. 

If one forms, namely, from (31) the linear combination 


(37) 


v- cJs&L = ir- a A Pic (p ^ ) 1 

obi U D (p) ° D' (Pfi) P - 


(38) 


the numerator of the left side of (38) represents a determinant which 
one obtains from the determinant of the characteristic equation (26) 
by replacing the k— th column by the C u ; thus one obtains a polynomial 

of at most m(n — l)— th degree in p. 


0=1 





(39) 


For (38) one can therefore write abbreviatedly 


n A (p) A '(p) pPr-1 A (Pn) 1 

Sfp) d(p) " ^o" d '(p^ P “ Pti 

while (37) obtains the form 

A|t (0) gfc* A„(Pn) Pu t 
Kinil - + f=Q P^DUP^) 


(38a) 


(37a) 


If one introduces, -moreover, for the homogeneous component solution according 
to (30) 


^Khom/^M^ 


mn-1 


-z: 

A=0 



and defines a function of time dependent only on the roots of the 

characteristic equation D(p) = 0 
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am— 1 

qp(t) 


3 p^ t 


k I=o~ Pn D ’(P|i) 

its A— th derivative with respect to time is 

.(A), 


<P 


(t) =2 


m=± p ^- 1 




g=?r D ’(Pn) 


( 4 o) 


(IjOa) 


and the solution of linear DES with constant disturbance functions and 
arbitrary starting conditions obtains the very clearly arranged form 


T A * (0> 

\ ‘ mr 


m(n-l) 

A=0 




mrt-1 

~K=Q 


(x) c K <p ft+l) (t) 


(in) 


whereas the derivative with respect to time of the K-th coordinate is 
given by the easily obtainable egression 


m(nr-l) 

A=0 




mn-1 

+izl 

A=0 


(x) C K cp (A+2) (t) 


( 4 la) 


S umm arizing, it may be repeated that the solution for constant 
disturbance terms may be represented in a clear form: the determinants 
D(p) are to be determined according to ( 26 ) and A K (p) with the 

coefficients ^b R according to the rule (39) whereas the coefficients 
and the function of time cp(t) are given by (30) and (10), 
respectively. 


3. Solution for Harmonic Disturbance Functions and Arbitrary 
Starting Conditions 


If, aside from the starting phenomena, periodic permanent conditions 
also are of Interest, It will be useful to write the disturbance functions 
in the form 


I„(t) = C^e 1 ^ 

may be a real as well as a caaplex quantity, 
the relation 


(36) yields as solution 
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x r ^ i ^ D " h 

K u=I 1 ^ 5 ~ ^ D ’ (?n) ift u - Pji tiar D * 


Kham/ P ^_ppt 


(Pji) 


(42) 


and one recognizes that particular - attention must he paid to cases where 


= P|1 


(43) 


Since all roots p^ were assumed different from, each other, this 
case will occur precisely for 

ifl = P 
x> *0 

If one writes (42) in the form 



-Po* 

0 O 

“ Po 


mn-1 A (Pu) 

fer D ’(pn) 


jiV- eV 


P H 


mn-1 D , (p,,) 

+ *S Khcm^' gu.t 

£=T !>•<*) 


(42a) 


and if now — > p Q , the first of the sum of (42a) assumes the indefinite 

form. 0/0, whereas the rest remains certainly finite when the real 
parts E(Pjj.) of the roots are less than zero. 

Thus the further deliberations concern only the first expression for 
which one has to perform the limiting process 


lim 

ifyj^ P Q 


iVL e Po t 
Xb - Po 



if one performs once more a limiting process with t 


00 
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r 


l±m tePot _ 
t — > ea 


0 for R(p o )<0 

os for R(p 0 ) > 0 


one recognizes that the solutions remain finite only if all real parts 
are R(pjj,)<0. 

Only under this assumption it seems logical to speak, in connection 
with (42), of a permanent condition for which then the relation 


n mn-1 

■^KDau** = ^ XZ 

. 0=1 p .=0 



ifip — p^ 


(46) 


is valid. 


From (31) follows for p = 1% the relation 


a fi.0,,1 mn— 1 A (Pu) 

^ ' __ ^ — UK 1 

p=0 (Py.) ^^3 — P[i 


and one obtains 


XicDau = e C ^D(i^')' 


( 47 ) 


b=l 


With the designations 




(47a) 


and 


D(ifl^) =ct^(fjy)+ iPuffiu) 


(47b) 


**Translator , s Note; The subscript "Dau" is an abbreviation of 
the German word "Eauerzustand" signifying "permanent." 
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respectively, and the fact that real as well as Imaginary part of (47) 
represent solutions, one obtains by transformation 


L KDau 


a u7'oit + 


0=1 


a •*> + 


COS — 


“"V ~ V 1 ?* Bin at 
“o * Pn 


1 


(48a) 


and 


L KDau 


y~ | a o7oK + sln + ^£ok — h£m CO s ^t 

0=1 j a^ 2 + p u 2 a^ 2 + f^ 2 

- 


(48b) 


respectively; however, for instance for (48a) one may write abbreviatedly 


‘•KDau “ ^ok cos (Got + e nK^ 
o=l 


(49) 


with the amplitude function 


V =, 

OK 


OK OK 

0^ 2 + P T 2 


(50) 


and the phase displacement 


= arctan 


/ q O S QK 

la r 

\0 OK 


__£o^ok\ 

+ P o 5 o J 


(51) 


One has to consider that the expression takes into account the excitation 
of each individual coordinate with a different frequency. One then has 
to stress not only that it is possible to eliminate the roots of the 
characteristic equation but also that one may successfully Include with the 
functions defined in ( 50 ) and ( 51 ) the influence of the o — th excitation 
on the K— th coordinate, and that one may obtain the total solution of 
the pe rmane nt equation by superposition of the partial solutions. 

With these statements we shall regard the most general case of linear 
DES with constant coefficients and n degrees of freedom as closed. 



30 


NACA EM 1197 


71. SOLUTION OF THE EQUATIONS OF MOTION 
DERIVED IN (IV) 


The application of this method to the treatment of the prescribed 
differential equation system ( 17 a), with the coefficients ( 17 b) for the 
disturbance functions 


and 


AS(t) = Const = ASmav 


AS(t) = ASmaxe 1 ^ 


proceeds particularly simply. 

Although the disturbance function ASmw v = const, is not exactly 

realizable in practice, it is well to note itB solution. It represents 
the upper limit for the controllability of an airplane by tabs and may 
serve as a comparative measure for control phenomena depending on time. 

The explicit representation for periodic disturbance functions may 
be particularly valuable when 1 nf ormation on resonance positions, amplitude 
functions, and phase displacements is required. 

If one wants, on the other hand, to estimate the possible error which 
may appear with the problem discussed in Section V,2 — to replace dis- 
turbance functions dependent on time by suitable selection of constant 
disturbance functions — it is suitable to use disturbance functions 

linear ily dependent on time. 

1. Solution for AS( T ) = Const. = AS^x 


If one takes into consideration that by the selection of the steady 
initial state all starting conditions are Aa=Arj = a = Tt = 0 and 
divides the equations of ( 17 &) by ASmmr. using further on the designations 


Aa _ y 
ASjuax ^ 
and m = 2 


and — = 3£i* one obtains according to (4l) with 
the relations 


n = 2 



* 


* 


A^O) 


'« D(O) 


2 


( A)-k m(A)^+^ 


\S*-f 


and 


X 


K 




(52a) 


respectively. 

far the determinant of the characteristic equation (26) one has to put with (23) 


D(P) - 


1 

( (0> % + (1) %P + S 2 ) ( (0)a ta + (1)a loP + (2)a taP 2 ) 
+0+o) ( (0, V + (1,a 2aP + P S ) 



arranged according to powers of p, the relations 


(53) 


U) 


H 
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(°>a 

f°)a 

la 


(°)a 

lil 

la 


to). 

1-q la 

a 0 - 

(0)a 

2tj 

(0) v 

& 1 = 

<°> S 

2ri 

H 

$ 

+ 

A 

o’ 

0 



/a\ /rt\ 

ke S, 


/-. \ /l \ 

Itj la 


1 (0 K 

fl-O 3S 

c. 

(0)a 2ri 1 

+ 

& 

H 

O 

+ 

o' 

p" 

s? 



(©a, <*>. 


1 


l(2) a 


1t| la 


la 


la 

a 3 = 

0 1 

+ 

0 <1)a 2a 

& 1 * 

0 1 


ar© valid for the coefficients of the determinant (53). 


J 


The determinant A K (p) vhich one obtains according to the rule (39) by replacing the rc-th 
column of (53) l>y "th© disturbance terms also arranging them, according to powers of p, is 




L 6 TZ W*I VOVH 





r, n /_\ 


^ 

2 =Z ”'' u kK - 

U=1 A=0 


with the coefficients 


(0)1,, 


( 0 ) 

C 1 “La. 


(°^u_ 


A, (0) 


( 0 ) 


"b2 


( 0 ), 


ft ln C 1 


(0), 


Co 


Ag(0) 


(1), 


( 1 ) 

C 1 a Lx 


C (l) a 

2a, 


(1), 


( 1 ) 


% C 1 


0 Or 


> 



r (2) a 
C 1 “let 


1 c x 

/n\ 

\ c --'h 1 =. 

n n 

^2 

(o\ m 
D2 = 

n 

- — 


f Rk\ 


\S~*f 


(54a) 


i 


u> 
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With (4o) the function cp(t ) 

3 

<?(t) = 3 

kL=0 


is represented, hy the relation 

eP^ T 

P^D’(Pji) 


(55) 


wherein are the simple roots of the characteristic equation (53) 

and with (33) 

D’(Pm) = “IT ^Up|jl-P„) (55a) 

olO 

For the numerical utilization of (55) one has to distinguish three 
cases . 

a) All roots real 

b) Two real and one pair con jugate-camp lex roots 

c) Two pairs conjugate-complex roots 

a) All roots of the characteristic equation real 

If all roots p^ are-real, the X— th derivatives of (55) appearing 
in (52) and (52a) are 


cp (X) (r) -ZI 

H=0 

and the solution (52) obtains immediately the form. 


Pn^Pt* 

D*(Pji) 


= 


Ak(0> 

D(0) 


M> 


V 


PnT 


(56) 


(57) 
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with 


and 


du - ± 


>,=0 


tj'Cph) 


x* = 21 

M-=0 


(57a) 


(57b) 


respectively. 


b) Two real and pne pair conjugate-ccanplex roots 

Since here and also in the case c) in the expression for cp^Or) 
complex quantities p^ appear, the representation of the solutions (52) 

in terms of real quantities proceeds somewhat Ibbs simply. However, the 
presence of conjugate-complex roots signifies only that the solution 
contains periodic functions. 


If one, therefore, splits off the two real roots p Q , p L and designates 

the pair of conjugate-complex roots ty p = k ±iv , one may write for 

2,3 1 1 

( 56 ) the expression 

cp^(r) = ^B o e' P 0 T + ^ c 0 e 5 l T + , (^b 1 sin Vjt + cos w- L r|e K l T (58) 

with 


<P ( 0 ) (t) = cp(r) and (0) B Q = B c 


etc. 
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and 


Bq = 


■ Po(p 0 - Pi) |(Po - Kl) 2 + n 2 ] 

1 

Pl(Pl - P^ j(Pl - *l) a + Vi 2 ] 

1 = ij* 1 ( P ° “ K l) ( p l ” *l) “ W l 2 ] + v l 2 ( p o + Pi - 2 ki) 

1 " “i*4 P ° “ “l) ( p l “ “l) - v l] - K l v l(Po + Pl “ 2«l)| 

(’‘l 2 + V l 2 ) [(fo - «l) 2 + V l 2 ] [f p l " + w l 2 ] 


N = w. 


(X) B q = p 0 ^B 


(X) , 

c o - Pl X °o 


(^) B _ (X~l)-D (X— 1) 

B i - "x B i _ *1 °1 


°1 “ V 1 B 1 + K 1 C 1 


Xf one now substitutes (58) into (52) and adds similar components 
of the solution, one obtains (58) in terms of real quantities 


X ic = d( 0) + a oK, eP ° + Pok® 5 - 1 - 1 * + v l T + PirCos v^e*! 7 ^ 


( 59 ) 
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with 


_ 2 _ 

>_ 

1=0 


a - Y~ ^h ^B 
a °K ~ 4-^ °K 


«K - f- ' X V X)b 1 


1=0 


> (59a) 


°H 


(X)v (>0 C 

d k o 0 


Sl*-i (X) Kk (X) 0i 


1=0 


and "by differentiation with respect to time the velocity 
‘ *b = ^oB® P ° T + PoK ePlT + |a lK sin w x t + FirCos w-jjje* 1 ' 1 ’ 


with 


a OK = p o <x ok 
^Ok = P 1^0K 


a lK = IC l a lK _ V lPlK 
P 1b = w l a lK + *l P lR 




c) Two pairs of conjugate-complex roots 

This case is solved, if one puts P c 1 = i^jiiv 0 
and p 2j 3 = K^ti-Wp with the relation 


(59b) 


(59c) 




U=0 


^B,*sin w_t + ^C„cos w, 


x> 


"U 1 


,i>e 


-K 0 T 


(60) 
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vith 



If one introduces here also (60) into ( 52 ) and forms the ’sum of similar 
terms, there appears here as the type of solution the form 


X, 


Ak(0) 

D(0) 


1 r 


+ 

VfcO 


K n T 


a UK sin w^t + (3 UK cos u 


vith 


a OK “ 


1=0 




p =3* ( x )h <*>C 

1=0 


(61) 


(6la) 
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and for the derivative 


*b = 21 <°x>K B±a w n T + p tK COB V 

P=0 1 


with 


“OK “ “ W oPuK 


” W t%K + R n* 3 UK 


2 . Solution for Harmonic Disturbance Function 

The solution of this special case is suitable for demonstrating the 
behavior of airplane controls by control surfaces and tabs for periodic 
excitations as they may be caused, for instance, by engine vibrations, 
fuselage-abending oscillations, or torsional oscillations of the stabilizing 
surfaces. However, one will have to pay special attention to the permanent 
forced oscillating state in comparison with the starting conditions. 

Noting the fact that the function e^^ 1 * which is dependent on time 
is independent of any summat ion, one obtains with (VjO and the relation 
(514-) as solution the permanent equation 




( 6 lc) 


j 


AS(t) = ZiSju a X e 


ittr 


e K b T 


( 6 lb) 


_ifiTA K (m) 
KDau " 0 D(m) 


(62) 


the values of the dete rmina nts A„(p) and D(p) have to be determined 

c •• 

according to ( 5 ^) and (53) for purely imaginary arguments p = ifi. 

Consequently it is useful to choose for these complex expressions 
the form 
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= 7 K (ft) + l& K (fl) 
D(ifl) = £(&) + it(fi) 




with the real functions of- fl 


y , (°>b _ ( 2 >b (l 2 

K K K 


s K = (1, b K n 


£ = a D - a 2 “ 2 - ^ 


\jr = a fl — a 

T 1 3 


(63) 


the coefficients of which are given with (53 a ) 3114 (5^)* respectively. 

If the real part of (6s) is to he regarded as solution, there results 
because of (38) 

X KDau = V 08 * 1 * + «r) W 


with the amplitude function 



V* * S K g 

i 2 + 1 2 


(6fe) 


and the phase displacement 


6 k = arc tan 1 - — — - ^ 

+ 

v. Hr ^ 


(64b) 
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3 . Solution for Linearily Dependent Disturbance Function AS(t) = At 


As already mentioned at the outset, this solution may serve for 
comparison of instantaneous and time— dependent sudden deflection of the 
tabs, if one superposes the two phenomena 


and 


AS(t) = ar 

for 

0 5 t < 0 = 

AS(r) = •< 

w 

0 

a(T- ti) 

for 

for 

0 ^ T < T. 
T 1 % T < c 


(65) 


As resultant deflection phenomenon one obtains for instance the disturbance 
function represented in figure 6 . Particular interest will he attached to 
that deflection phenomenon which characterizes the maximum velocity of 
deflection of the tabs possible in practice and is given by 


a = ^ max 
T 1 


( 65 a) 


To attain the solution one will go back to the general form of eguatlon 
(36). 

With the fact already stressed in Section VI, 1 that all starting 
conditions are zero, there follows with 


»„Ct> - 

T 1 

after performance of the integration for the ratio -^ 2 — = x-j_ 

ASmax 

and -^3 — = Xq; respectively, 

ASjnax 


v _ 1 

“ 77 2 d a 

'A u-Ti „ A 


|e^ T 



( 66 ) . 
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With the relation (64), however, one can write 



|a*(o) 

1 

[D(0) 


± A >) 

h= 0 ph 2d ’(ph) ^OP^’CPtx) 


(66a) 


or, after according rearrangement and introduction of the notation in 
powers for 


1 

< 

^ (0) T + 4 - (ok 

e^ T 1 

T 1 

0 ( 0 ) ■ $0 

Pn 2 P' (P^) P^' (Pjj.) 


3 (2 , b _ 3 S 

^=0 K E'(P(i) ^i=0 X=1 ' 

If one farms from (55) the integral 


4- (U„ e p)lT 


H=0 P^P^Pn) 


>(66b) 


r 3 

/ cp(£)dl = Y~ ■< 
1 0 

e p t J,T 

> 

1 

P^P’(Pn) 

Pp^'Cp^) f 


L J 


m 


and introduces the designation / cp(E )di = cp^ ^(t) one may write for 

Jo 

(66), in analogy to (52), also 


k - i +± W V (W »(r) -£ £ 4^— i 

1 |P(0) \=0 n=0 X=1 p ti 2 D»(p H ) I 


(66c) 


Corresponding to the three combinations of the roots of the characteristic 
equation (53) treated in Section VI, 1 these three possibilities could be 
investigated here also. 
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^3 


Let It suffice to point out the procedure for the case of two pairs 
of conjugate— complex roots. As analogon to (60) one may select the 
expression 



\ 

’ 

f 9 ( 1)46 = 2Ei 1 

[ (_i) B u sin w^t + (_l) C D cos w uT ] e K P T - (-l) C u 

- 

/ n=0 

'0 


with 


(67a) 


(-D 


B o " 


“T 


'u 


+ w u C u 
+ v 0 2 


(-D 


'v 


'SCp 

K 2 . v 2 


with the B-y and 0 o to he calculated according to (60h). With (67a) 
and (6o)j (66c) is transformed into the form 






AJO) 


T 1 


D(0) 


•T + X P 

i>=o L 


a sin w t + B cos w t 
UK O k UK X> 


K X> T 


( ^Co + ^ 1 ^Ci + 


2 (X)v 


H=0 1=1 p/Wp^J 


> 


(-58) 


with 



> (68a) 
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If one introduces into (68) according to rule (65) the argument (t - tj) 

instead of the independent variable t and superposes the solutions thus 
obtained, the solution for the deflection phenomenon sketched in figure 6 
assumes the shape 


Mr "N 


a k (0) 


+ > 


■^- T l D(0) 'tgo 


K„T 


a m sin w u t + 0 UK cos e ^ 


(69) 


with 


a UK = ~ ja-\)ic(l - e K n T lcos w^yr^ + p^e K u T lsin w^t-j^ 

Pw = ^| a UK e ~ Kt>11al11 w n T l + " e“ K o T l cos w^j 


> 


whereas for the derivative the relation 


i K - x{ a n K sln v n T + Kk cob v| eKuT 


1>=0 ^ 


with 


%K = — W, oPuiC 


= + K ul 5 UK 


> 


exists and one must note that (69) is valid only in the 

T-l < T < 00 . 


(69a) 


(69b) 


(69c) 


interval 



KACA 524 1197 


^5 


The first partial report may "be closed vith these theoretical exposi- 
tions. In the second partial report the results of the present report 
will he applied to large aircraft controlled by tabs. 


YU. SUMMARY 


The present report is intended as a contribution toward the clarification 
of problems arising for longitudinal motions of airplanes which are con- 
trolled not by direct control surface activation but indirectly by tabs. 

After a general derivation of the equations of motion for the four 
degrees of freedom X , Z , -8 , and q and taking into consideration that 

o o 

the control surfaces are not weight— compensated, permissible assumptions 
are introduced so that the system of four degrees of freedom is reduced to 
a system with the two degrees of freedom a and T| . By the 

substitution t = t \Jq a representation of the equations of motion which 
is independent of the dynamic pressure becomes possible Inasmuch as one 
regards the dynamic pressure ratio qg/q an d the downwash factor c&/cta, as 

constant . 

While stability and oscillation conditions may be discussed already 
with a knowledge of the solution of homogeneous differential equation 
systems, one must know also the solutions of inhomogeneous differential 
equation systems for the investigation of starting phenomena . Since the 
required mathematical expedients appear suitable not only for the solution 
of the prescribed problem but also for the analytical description of 
similar problems, they are treated in the most general formulation for a 
linear differential equation system of n degrees of freedom with constant 
coefficients, and arbitrary starting conditions and disturbance functions 
as far as the latter are integrable. 

From the general solution special solutions for constant and periodic 
disturbance f unctions are developed. 

With the definition of a function cp(t) which depends only on the 
roots of the characteristic equation one may obtain for constant disturbance 
functions a solution of a very clearly arranged form; it consists of the 
three constituents: constant, inhomogeneous, and homogeneous part of 

the solution. 

Far periodic excitations one determines first the conditions which 
are required to make a permanent state possible. Then one can demonstrate 
that amplitude and phase— displacement functions exist which include the 
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influence of the o— th excitation on the K—th coordinate. The total 
solution of the permanent equations is then obtained by superposition of 
the partial solutions. 

These results yield, applied to the prescribed problem, the desired 
solutions] with regard to numerical calculation real representation of tha 
functions is particularly emphasized. 


Translated by Mary L. Mahler 
National Advisory Committee 
for Aeronautics 
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Figure 2.- Direct tab control. 



Figure 3.- Indirect tab control. 
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